CONVERGENCE ANALYSIS OF A COLOCATED FINITE VOLUME 
SCHEME FOR THE INCOMPRESSIBLE NAVIER-STOKES 
EQUATIONS ON GENERAL 2 OR 3D MESHES 
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Abstract. We study a colocated cell centered finite volume method for the approximation 
of the incompressible Navier-Stokes equations posed on a 2D or 3D finite domain. The discrete 
unknowns are the components of the velocity and the pressures, all of them colocated at the center 
of the cells of a unique mesh; hence the need for a stabilization technique, which we choose of the 
Brczzi-Pitkaranta type. The scheme features two essential properties: the discrete gradient is the 
transposed of the divergence terms and the discrete trilinear form associated to nonlinear advective 
terms vanishes on discrete divergence free velocity fields. As a consequence, the scheme is proved 
to be unconditionally stable and convergent for the Stokes problem, the steady and the transient 
Navier-Stokes equations. In this latter case, for a given sequence of approximate solutions computed 
on meshes the size of which tends to zero, we prove, up to a subsequence, the L 2 -convergence of 
the components of the velocity, and, in the steady case, the weak L 2 -convergence of the pressure. 
The proof relics on the study of space and time translates of approximate solutions, which allows 
the application of Kolmogorov's theorem. The limit of this subsequence is then shown to be a weak 
solution of the Navier-Stokes equations. Numerical examples are performed to obtain numerical 
convergence rates in both the linear and the nonlinear case. 
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1. Introduction. We are interested in this paper in finding an approximation 
of the fields u = (u (l) ) 4= i,...,d : f2 x [0, T] — > R d , and p : fl x [0, T] -»■ R, weak solution 
to the incompressible Navier-Stokes equations which write: 



<9 t u w - vAu il) + d t p + ^2u U] ] d 3 u {l) = in Q x (0, T), for i = 1, . . . , d, 

(!•!) 

divu = J2 d ^ il) =0 in fix (0,T). 



nstocontt 



with a homogeneous Dirichlct boundary condition for u and the initial condition 
(1.2) 0) = in fl for i = 1, . . . , d. 

In the above equations, vy\ i = 1, . . . , d denote the components of the velocity of a 
fluid which flows in a domain fl during the time (0, T), p denotes the pressure, v > 
stands for the viscosity of the fluid. We make the following assumptions: 



(1.3) 
(1.4) 
(1.5) 
(1.6) 
(1.7) 



fl is a polygonal open bounded connected subset of R d , d 
T > is the finite duration of the flow, 
v £ (0, +oo), 
u ini € L 2 (Q) d , 



2 or 3, 



/« e L 2 {fl x (0,T)), for i = 1, 



,d. 
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We denote by x = (a;^)j=i ... d an y point of Q, by |.| the Euclidean norm in 

d 

M. d , i.e.: \x\ 2 = ^^(a;^) 2 and by da; the d-dimensional Lebesgue measure da; = 

i=l 

dx^-.-dx^. 

The weak sense that we consider for the Navier-Stokes equations is the following. 

Definition 1.1 (Weak solution for the transient Navier-Stokes equations). 
Under hypotheses ll.3\) - [L~7\) , let the function space E(£l) be defined by: 



(1.8) 



E(Ct) := {v = (w W ) i= i,...,d G fl£(fi) d ,divu = a.e. in fi}. 



Then u is called a weak solution of iTO]) - ^) ifu G L 2 (Q, T; £(ft))nL°°(0, T; L 2 (fl) d ) 



id: 



( Vtp G L 2 (0,T;E(n))C)C™(n x (-oo,T)) d , 

T 



(1.9)1 



o Jn 



l(x,t) ■ dt(p(x,t)dxdt — / Ui n i(x) ■ (p(x, 0) da; 



o Jn 



Vu(x,t) : Vip(x,t)dxdt + 



b(u(-,t),u(;t),^(-,t))dt 

T 

f(x) ■ <p{x, t) dx dt 



o Jn 



where, for all u,v G HQ(fl) d and for a.e. x G f2, we use the following notation: 



V«(x) : Vu(x) = ^Vu w (x) • Vu w (a 



eO 



nstocontf t 



and where the trilinear form b(., ., .) is defined, for all u,v,w G (HQ(£l)) d , by 



d d 



[1.10) 



b(u,v,w)=Y^T f u^(x)d i v^(x)w ( ~ k) (x)dx. 

1 1 „' 1 J £1 



fc=li=l 



def tricont 



Jnst. prnTit.ft. In St 

Remark 1.1. From U.9\) . we get that a weak solution u of ll.l\) - kLty) in the sense 
of Definition[Tj\ satisfies d t u G L 4 / d (0, T; E(Q,)'), and is therefore a weak solution in 
the classical sense, such that u(-,0) is the orthogonal L 2 -projection of Uj n ; on {v G 
L 2 (i}) d , divw = 0,trace(w ■ ngn,dfl) = 0} (see for example f3bj or J7$). 

Numerical schemes for the Stokes equations and the Navier-Stokes equations have 
been extensively studied: see [23l [33l [3U [35l [24] and references therein. Among 
different schemes, finite element schemes and finite volume schemes are frequently 
used for mathematical or engineering studies. An advantage of finite volume schemes 
is that the unknowns are approximated by piecewise constant functions: this makes it 
easy to take into account additional nonlinear phenomena or the coupling with alge- 
braic or differential equations, for instance in the case of reactive flows; in particular, 
one can find in [33j the presentation of the classical finite volume scheme on rect- 
angular meshes, which has been the basis of many industrial applications. However, 
the use of rectangular grids makes an important limitation to the type of domain 
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which can be gridded and more recently, finite volume schemes for the Navier-Stokes 
equations on triangular grids have been presented: see for example where the 
vorticity formulation is used and [5] where primal variables are used with a Chorin 
type projection method to ensure the divergence condition. Proofs of convergence for 
finite volume type schemes for the Stokes and steady-state Navier-Stokes equations 
are have recently been given for staggered grids [9], [26], [13], [H], [4], following the 
pioneering work of Nicolaidcs et al. [31] , [32] . 

In this paper, we propose the mathematical and numerical analysis of a discretiza- 
tion method which uses the primitive variables, that is the velocity and the pressure, 
both approximated by piecewise constant functions on the cells of a 2D or 3D mesh. 
We emphasize that the approximate velocities and pressures are colocated, and there- 
fore, no dual grid is needed. The only requirement on the mesh is a geometrical 
assumption needed for the consistency of the approximate diffusion flux (see |15| and 
section @ for a precise definition of the admissible discretizations). 

As far as we know, this work is a first proof of the convergence, of a finite volume 
scheme which is of large interest in industry. Indeed, industrial CFD codes (see e.g. 
[28| . [1]) use colocated cell centered finite volume schemes; leaving aside implemen- 
tation considerations, the principle of these schemes seems to differ from the present 
scheme only by the stabilization choice. The main reasons why this scheme is so 
popular in industry are: 

• a colocated arrangement of the unknowns, 

• a very cheap assembling step, (no numerical integration to perform) 

• an easy coupling with other systems of equations. 

The finite volume scheme studied here is based on three basic ingredients. First, 
a stabilization technique a la Brezzi-Pikaranta [8] is used to cope with the instability 
of colocated velocity /pressure approximation spaces. Second, the discretization of 
the pressure gradient in the momentum balance equation is performed to ensure, 
by construction, that it is the transpose of the divergence term of the continuity 
constraint. Finally, the contribution of the discrete nonlinear advection term to the 
kinetic energy balance vanishes for discrete divergence free velocity fields, as in the 
continuous case. These features appear to be essential in the proof of convergence. 

We are then able to prove the stability of the scheme and the convergence of dis- 
crete solutions towards a solution of the continuous problem when the size of the mesh 
tends to zero, for the steady linear case (generalized Stokes problem), the stationary 
and the transient Navier-Stokes equations, in 2D and 3D. Our results are valid for 
general meshes, do not require any assumption on the regularity of the continuous 
solution nor, in the nonlinear case, any small data condition. We emphasize that the 
convergence of the fully discrete (time and space) approximation is proven here, using 
an original estimate on the time translates, which yields, combined with a classical 
estimate on the space translates, a sufficient relative compactness property. 

An error analysis is performed in the steady linear case, under regularity assump- 
tions on the solution. An error bound of order 0.5 with respect to the step size is 
obtained in the discrete H 1 norm and the L 2 norm for respectively the velocity and 
the pressure. Of course, this is probably not a sharp estimate, as can be seen from 
the numerical results shown in Section [S] Indeed, a better rate of convergence can be 
proved under additional assumptions on the mesh 20\. 
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This paper is organized as follows. In section ® we introduce the discretization 
tools together with some discrete functional analysis tools. Section [3] is devoted to 
the linear steady problem (Stokes problem), for which the finite volume scheme is 
given and convergence analysis and error estimates are detailed. The complete finite 
volume scheme for the nonlinear case is presented in section [4] in both the steady and 
transient cases. We then develop the analysis of its convergence to a weak solution 
of the continuous problem. We give some numerical results in section and finally 
conclude with some remarks on open problems (section [5]) . 



2. Spatial discretization and discrete functional analysis. 

2.1. Admissible discretization of f2. We first recall the notion of admissible 
discretization for a finite volume method, which is given in |15j . 

Definition 2.1 (Admissible discretization, steady case). Let fl be an open 
bounded polygonal (polyhedral if d = 3) subset ofM. d , and dil = fi\fi its boundary. An 
admissible finite volume discretization ofQ, denoted by T>, is given by T> = (Ai, £, V), 



- Ai is a finite family of non empty open polygonal convex disjoint subsets of 
n (the "control volumes") such that Q = UkemK- For any K £ Ai, let 
dK = K\K be the boundary of K and tsik > denote the area of K. 

- £ is a finite family of disjoint subsets of £1 (the "edges" of the mesh), such 
that, for all a £ £, there exists a hyperplane E of K d and K £ Ai with 
a = dK n E and a is a non empty open subset of E. We then denote by 
m CT > the (d-1)- dimensional measure of a. We assume thatjor all K £ Ai, 
there exists a subset £k of £ such that dK = U a ££ K a. It then results from 
the previous hypotheses that, for all a £ £, either a C dfl or there exists 
(K, L) £ Ai 2 with K L such that K n L = ~o; we denote in the latter case 
a = K\L. 

- V is a family of points of f2 indexed by M., denoted by V = (xk)kgm- The 
coordinates of xk o,re denoted by xS , i = 1, . . . , d. The family V is such that, 
for all K £ Ai, xk €E K. Furthermore, for all a £ £ such that there exists 
(K, L) £ M 2 with a = K\L, it is assumed that the straight line (xk, xl) going 
through xk and xj, is orthogonal to K\L. For all K £ Ai and all a £ £k , let 
z G be the orthogonal projection of xk on a. We suppose that z a £ a. 

An example of two neighbouring control volumes K and L of M. is depicted in 



The following notations are used. The size of the discretization is defined by: 



For all K £ M. and a £ £k, we denote by n^- jCr the unit vector normal to a outward 
to K. We denote by dK,a the Euclidean distance between xk and a. The set of 
interior (resp. boundary) edges is denoted by & m t (resp. £ e xt), that is & m t = {cr £ £; 
a <f_ dfl} (resp. £ C xt = {o - £ £ ', o~ C dfl}). For all K £ M, we denote by Mk the 
subset of Ai of the neighbouring control volumes. For all K £ Ai and L £ Mk, we 
set thkl = n-K,K\L, we denote by d K \h the Euclidean distance between xk and xl- 
We shall measure the regularity of the mesh through the function regul(2?) defined 



secdisc 



where: 



size(£>) = sup{diam(A'), K £ Ai}. 



by 



(2-1) 




"I 



dK,K\L 
d-K\L 




card(£ k) 



1 



K £ Ai \ . 



i 



Fig. 2.1. Notations for an admissible mesh 
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2.2. Discrete functional properties. Finite volume schemes are discrete bal- 
ance equations with an adequate approximation of the fluxes, see e.g. |15) . Recent 
works dealing with cell centered finite volume methods for elliptic problems [21] , [16] , 
|14j introduce an equivalent variational formulation in adequate functional spaces. 
Here we shall follow this latter path, also introducing discrete analogues of the con- 
tinuous Laplace, gradient, divergence and transport operators, each of them featuring 
properties similar to their continuous counterparts. 

Definition 2.2. Let be an open bounded polygonal subset ofW. d , with deN,. 
Let T> = (Ai,£,V) be an admissible finite volume discretization of Q in the sense 
of definition \2.1l We denote by H-r>{Q) C L 2 (fl) the space of functions which are 
piecewise constant on each control volume K G M.. For all w G Hx>{Q) and for 
all K G M, we denote by u>k the constant value of w in K. The space Hx>(il) is 
embedded with the following Euclidean structure: For (v,w) G (Hd(CI)) , we first 
define the following inner product ( corresponding to Neumann boundary conditions ) 



oo 



(2.2) 



(v, w)v 



- E 

2 ^ 

KeM LeNf 



E 



™k\l 
dx\L 



{v L - v K )(w L - w K )- 



amundis 



We then define another inner product ( corresponding to Dirichlet boundary conditions ) 

m„ 



(2.3) 



[v,w\v = {v,w)j> 



E E 



-v K w K - 



lK,cr 



Next, we define a seminorm and a norm in Ht>{^1) (thanks to the discrete Poincare 
inequality \2.4\ ) given below) by 

\w\v = ((w,w) v ) 1/2 , \\w\\ v = {[w,w] v ) 1/2 . 

We define the interpolation operator Pp ■ C(f2) — ► Ht>(Q) by (Pd^>)k = <p(xk), for 
all KeM, for all ip G C(fi). 

Similarly, for u = (u«) i= i,... )d G {H v (n)) d , v = (u«) i= i,..., d G (H v {Q)) d and 
w = (iy w ) J= i,...,d G (H- D (fl)) d , we define: 



amudis 




\V,W\v 



weakconvgrad 



andP v : C(n) d ^ H v (n) d by (P v ip) K = w(x^), for all K e M, for all (p e C(n)° 
The discrete Poincare inequalities (see |15j ) write: 

(2.4) IMU 2 (o) < diam(n)||iu||i,, Vu; G H V {Q), 
and there exists Cn > 0, only depending on f2, such that 

(2.5) IIHli*(n) - c o\w\t» Vw e #z>(n) with ^ w(x)da; = 0. 

We define a discrete divergence operator div-p : (ifp(Sl)) d — > ifp(O), by: 

(2.6) divp(u)(iE) = — ^— A A -L • (u K + u L ), for a.e. x £ K,WK e 7W, 



m A - 



with 
(2.7) 



Akl = ™ KlL XL „ = J m K|i iikx, WT eMVLe Mr- 



We then set Eb(fi) = {«£ (H v (ty) d , divp(w) = 0}. 

Remark 2.1. ^4ny definition of A^l such that A^l = m i<"|L^ifL n ArL with okl > 
tm<i okl + ilk = 1, combined with the definition <Hv-d(u)(x) = ^LeM K ^ KL ' 
uk — Alk • ul), produces the same results of convergence as those which are proven 
in this paper. On particular meshes, one can prove a better error estimate, choosing 
okl = d(xi,, K\L) / duL (see \20i ). Nevertheless, in the general framework of this 
paper, other choices do not improve the convergence result and the error estimate. 
Therefore, we set in this paper okl = 1/2, which corresponds to \2. 7\ ). The advantage 
of this choice is that it leads to simpler notations and shorter equations. 

The adjoint of this discrete divergence defines a discrete gradient Vx> : -ffx>(f2) — ► 



(2.8) 



(Vt>u)k = — V A KL {u L - u K ), VK eM, Vue H V {Q). 

TY1 Ty ' * 



This operator Vx> then satisfies the following property. 

Proposition 2.3. Let (2?("^) m eN be a sequence of admissible discretizations of 
in the sense of Definition [2~H such that lim sizc(T> i ' m ^) = 0. Let us assume that 

m — »oo 

there exists C > and a G [0, 2) and a sequence (u''"') me n such that vy 71 ' G H V ( m ) (fi) 
and |u (m) |f> m < C size(V^)- a , for all m G N. 
Then the following property holds: 

(2.9) lim / (Pp <p{x)V v u^ n] (x) + u^ n] (x)V ip(x)) dx = 0, Vtp G C c °°(fi), 

rn->+caJ n V / 

and therefore: 



(2.10) 



lim 

m — >+oo 



U ^{x) ■ p Vm ii){x)Ax = o, w g c c °°(fi) d n E(n), 



poindis 



poindismoy 



divdisc 



defAB 



discgrad 



wevgrad 



wevdiv 



where E(Q) is defined by U.8\) . 



cpct 



Proof. Let us assume the hypotheses of the above lemma, and let i = 1, . . . , d and 
tp G C£°(n) be given. Let us study, for m G N, the term 



£grad 

|), we get that 



MafAB Miscgrad 

From dUTJ) and (gU), we get that 



where 



r( kl = ( o(v( x k) + <p{xl)) — [ <p(x)drf(x) ) n KL . 

\ 2 m K\LJK\L J 



Thanks to the Cauchy-Schwarz inequality, 

One has X^o-ef- t a=K\L m K\LdKL < cZm(Q). Thanks to the existence of C v > which 

only depends on tp such that < C^sizc and since a < 2, we then get 

that 

lim T,4 m) = 0, 

which yields (|2J])nj 

PROPOSITION 2.4 (Discrete Rellich theorem) . Let (P'"^) mG N Z>e a sequence of ad- 
missible discretizations oj m i/ie sense of definition \2.1\ such that lim size(2?' m - 1 ) = 

m — >oo 

0. Lei us assume that there exists C > and a sequence (u m )wien sweft £/iai 
u (m) G % m) (0) and ||w (m) ||p m < C 1 /or a?/ m G N. 

Then, there exists u G i?o(f2) anc ^ a subsequence of («' m ')meN, again denoted 

J. i/ie sequence («' m ') mE N converges in L 2 (f2) to u as m — > +oo, 
£. /or a/Z 93 G C£°(fi), we Ziawe 

(2.11) lim [u {m \P v ip] v = I VS(x) • V^(x)dcc, 

m->+oo 7 n 

hjr. wT-ad 

5. Vp m u < -'™' weakly converges to Vu in L 2 (f2) as m — > +oo and l£. 9\) holds. 

Proof. The proof of the first two items is given in [15] (see proof of Theorem 
91. pp 773-774). Since we have \u (m) \v m < ||« lm; ||x> m ) we can apply proposition 12.31 
which gives the third item. □ 

Remark 2.2. Following fWj , if we denote 

T> K ,a = {tx K + (1 - t)y, t G (0, 1), y G a}, MK G M., Vcr G £ K , 
we maw alternatively define a discrete gradient Vx> : Hj)(fl) — > (L 2 (f2)) d , ow: 
/or aZZ A'gM, 

Vp«(a;) = d7tj;( UL ~ u k) u kl, for a.e. x G V K , K \ L U V L ^ K \ Ll VL G AZ*r, 
Vpu(i) = -j (0 - UK)riK.a, for a.e. x G Pr>, Vcr G £sr n £ e xt- 

<*K> 



cstd 



secf vslin 



A result similar to that of Proposition \z^\ holds with this definition of a discrete gradi- 
ent, and in fact, it can be shown that the weak convergence o/V^u'" 1 ' is equivalent 
to the weak convergence of Vx) m u^ m ^ . 

3. Approximation of the linear steady problem. 

3.1. The Stokes problem. We first study the following linear steady problem: 
find an approximation of u and p, weak solution to the generalized Stokes equations, 
which write: 



(3.1) 



rju — i/Au + Vp = / in 
divu = in Q, 



stocont 



For this problem, the following assumptions are made: 



(3.2) O is a polygonal open bounded connected subset of M , d — 2 or 3 



hypomega 



(3.3) 



v e (0,+oo), 77 e [0,+oo), 



hypnu 



weaksol 



(3.4) 



We then consider the following weak sense for problem 

Definition 3.1 (Weak solution for the steady Stokes equations). 

Under hypotheses &3.2\) - (3~4\ l, let E(Q) be defined by il.8\) . Then (u,p) is called 

a weak solution of \3. 1\) (see e.g. \36f or [7|/) if 



hypfg 



(3.5) 



u e E(Q), p G L 2 (fl) with J n p(x)dx = 0, 

77 / u(x) ■ v(x)dx + v I Vu(x) : Wv(x)dx— 
n Jn 

p(x)divv(x)dx = I f(x)-v{x)dx, VveH^{n) d . 



stocontf 



The existence and uniqueness of the weak solution of ( |3.ip in the sense of the 
above definition is a classical result (again, see e.g. or [7]). 

3.2. The finite volume scheme. Under hypotheses (|3.2|l - (|3.4p . let T> be an 

admissible discretization of Q in the sense of Definition 12. ll It is then natural to write 
an approximate problem to the Stokes problem (|3 . 5|) in the following way. 



ueBp(fl), p e H v (n) with / p(x)d 



(3.6) 



x = 



rj / u(x) ■ v(x)dx + v\u, v]x> 

- / p(x)div v (v)(x)dx = I f(x)-v(x)dx VveH v (n) c 
Jn Jn 



schvf 



As we use a colocated approximation for the velocity and the pressure fields, the 
scheme must be stabilized. Using a non-consistant stabilization a la Brczzi-Pitkaranta 



seccvglin 



estlll2 



, we then look for (u,p) such that 

(u,p) G H v (n) d x H V (Q) with f p(x)dx = 

Jn 



(3.7) { 



T) I u(x) ■ v(x)dx + v\u, v]d 
p(x)divx>(v)(x)dx 



f(x) ■ v(x)dx 
divx>{u)(x)q(x)dx = —A size(V) a (p,q)j> 



W 6 H v (n) d 

v q e H v (n) 



where A > and a G (0,2) are adjustable parameters of the scheme which will have 
to be tuned in order to make a balance between accuracy and stability. 

System (|3.7|l is equivalent to finding the family of vectors (uk)k£M C R d , and 
scalars (pk)kgm C K solution of the system of equations obtained by writing for 
each control volume K oi Ai: 



T) Vl\k UK ~ V 



E m K\L , . \ - 



(0-u K ) 



(3.8) J 



V A KL {pl-Pk)= \ f(x)da 

TZr.. JK 



LeM K 



A KL -(u K +u L )-Xsize(V) a ^j^-(Pl-Pk)=0 



supplemented by the relation 
(3.9) 



^ m A ' PA" = 



KEM 



Defining p a = (pk + Pl)/% if c = K\L, and p CT = px if <r G £ C xt H £a, and using the 
fact that X)o-ef: A - m cr n A",o- = 0, one notices that: J2l£AT k Akl (pl — Pk) is in fact 
equal to X^o-gf ^<yP<y^K,<y, thus yielding a conservative form, which shows that (|3.8I) 
is indeed a finite volume scheme. 

The existence of a solution to (|3.7p will be proven below. 

3.3. Study of the scheme in the linear case. We first prove a stability 
estimate for the velocity. 

Proposition 3.2 (Discrete H 1 estimate on velocities). Under hypotheses 
\3.4\ ), let T> be an admissible discretization of in the sense of definition \2.1\ Let 
A G (0, +oo) and a G (0,2) be given. Let (u,p) G Hx>{Q) d X Ht>{£1) be a solution to 
\3. 7| ). TTien the following inequalities hold: 



schvf 



schvf S 



moypnulle 



fS 







(3.10) 

and 

(3.11) 



v\\u\\ v < diam(f7)j|/|| (i2(n)) d, 
v A size(2?) Q \p\ 2 v < diam(r>) 2 ||/||2 i2(0))d . 



estimu 



est imp 
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Proof. We apply (|3.7p setting v = u. We get 



j] I u{x) 2 dx + Hl u llx> 



p(x)divx>(«)(x)da; = / /(x) ■ v(x)dx 



Since r/ > 0, the second equation of (|3.7[) with q = p and Young's inequality yield 
that: 

/ u(x) 2 dx + ^||w||f, + Asizc(:D) a |p|f, < 

dlam (") II f II 2 | " |L,||2 

2l/ H/II(LW + 2diam(0)2 ||U|l (i 2 (n))-*- 



Using the Poincare inequality (|2.4j) gives 



exunpos 



estp 



i>Hlf>+A sizc(P) Q IpH, < 
which leads to (f3TT0|) and (|3TTT|) . □ 



diam(fi) 2 
2^ 



l/ll 



(L 2 (n)) d + 2" u 



We can now state the existence and the uniqueness of a discrete solution to (|3.7p . 

Corollary 3.3. [Existence and uniqueness of a solution to the finite 
volume scheme] Under hypotheses iS.Zfy - pTfy , letT> be an admissible discretization 
of fi in the sense of Definition [KT[ Let X £ (0, +oo) and a £ (0,2) be given. Then 
there exists a unique solution to \3. 7\ ). 

Proof. System (|3.7I) is a linear system. Assume that / = 0. From propositions 
3.21 and using (|2.5p . we get that u = and p = 0. This proves that the linear system 
(|3.7j) is invertible. □ 

We then prove the following strong estimate on the pressures. 

PROPOSITION 3.4 (L 2 estimate on pressures). Under hypotheses H3.ty) - (3~J\ ), let 
T> be an admissible discretization of Q in the sense of definition \2.1\ and let 9 > 
be such that rcgul(P) > 9. Let A £ (0, +oo) and a £ (0,2) be given. Let (u,p) £ 
i?x>(0) x iZ-p(f2) be a solution to \3. 7\ ). Then there exists Che,sQV>ly depending on d, 
£1, r\, v , X, a and 9, and not on size(P), such that the following inequality holds: 



(3.12) 



\p\\l*{si) < <^\\f\\(L^n)y- 



estimsp 



Proof. We first apply a result by Necas [55]: thanks to J n p(x)dx = 0, there exists 
C 2dei'rUiT. which only depends on d and f2, and v £ Hq (fl) d such that divv(x) = p(x) 



for a.e. x £ f2 and 



(3-13) IMIffi(n)<i < Cj2]IMU 2 (n)- 

We then set 

= — [ 5 (l) (x)d 7 (x), Va££, Vi = l,...,d. 

m CT J a 

(note that vy = for all a £ £ CKt and i = 1, . . . , d) and we define v £ Hx>(Q) d by 



W 

,j k = 



— [ v {i) (x)dx, VK £ M, V* = 1,. ..,d. 
mx J K 



' K 
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Applying the results given p 777 in [T5], we get that there exists C jjdi^nli) , only 
depending on d and (9, such that 



(3.14) 

and 
(3.15) 

We then have 



(t^-tW^Cfe, 



diam(A) 



(Vu W (x)) 2 dx, 



1Mb < c^\v\\ H i [n)d . 

I p(x)div v v(x)dx = V] Pa: V] ^4ifL • («x + «l) = Tbdrp^fW 



kem LeM K 



where 



= £ PK X] ' 

KeM LeM K 



^2 PK / v ^ ' n ^i d 7(a;) 

A A J r- J\[jf K\L 



KeM he 



p(x)divv(x)dx = \\p\\ 2 2(n) , 



and 



KeM LeN K 

= E 

a=K\Le£i, 

We then have, thanks to the Cauchy-Schwarz inequality 



T [3] = E PK E mK \ L (\( v k + v l) -v k \ l J -n KL 
M LeM K ^ ' 

E m A"|L(PK - Pl) I ^("K + v L ) - lf K | L 1 ■ n KL . 



\p\v E m K|L d ifL ( \{ V K + V L ) - V K]L 



cr=K\Le£ in 



inebookl 



continjhld 



Applying Inequality (|3. 14[) and thanks to (|(% + Vl) — vk\l) 2 < \{{ v k — vk\lY 
{vl - v k \l) 2 ), we get that 



%<\p\l £ d KL Cfosize(P) [ Y,(Vv^(x)) 2 dx 



rr=K\LeS iB 

This in turn implies the existence of CU e ^t^pOnly depending on d and 6, such that 

T|<q4]si Z e(I?) 2 b||,^||| ol(n)li . 
Thanks to (|3.13[) . we then get, gathering the previous results 



(3.16) 



p{x)div v v(x)dx > |b||| 2(n) - Cfaawe(V)\ P \ v Cft\\p\\ L 2 {{l) . 

n 



latche 



We then introduce v as a test function in (|3.7p . We get 



cvgce 



(3.17) / p(x)divx>(v)(x)dx = rj j u(x) ■ v(x)dx + u[u, v]x> — / f(x)-v(x)dx. 
Jn Jn Jn 

Applying the discrete Poincare inequality, (|3. 15[1 and (|3.16|) . we get the existence of 



vdansschvf 



C 5 QS ^5^3iepending on d, f2, /, r), u, A and 6, such that 

\\p\\h(n) - q5pize(D)|p|2><^||p|| L 2 (n) < qg (U\t> + ll/IUw) INU»(n). 

We now apply (f37TO|) and (f37TT]) . Since sizc(P) 2 < sizc(P) Q diam(fi) 2 - Q , the condi- 
tion a < 2 suffices to produce (|3.12p from the above inequality, a factor 1/A being- 
introduced in the expression of Cjj] (it is therefore not possible to let A tend to in 

Wc then have the following result, which states the convergence of the scheme 

Proposition 3.5 (Convergence in the linear case). Under hypotheses L3.2\ 
PIP , let (u,p) be the unique weak solution of the Stokes problem i3.1]) in the sense 
of definition \3.1\ Let A G (0, +oo), a G (0,2) and 9 > be given and let T> be an 
admissible discretization of fl in the sense of definition \ 2. ?1 such that rcgul(2?) > 9. 
Let (u,p) G H v (n) d x H v (ty be the unique solution to \3. 7| ). 

Then u converges to u in (L 2 (Q)) d andp weakly converges to p in L 2 (Q) as size(2?) 
tends to 0. 

Proof. Under the hypotheses of the above proposition, let (2?( m )) m€ pj be a se- 
quence of admissible discretizations of Q in the sense of definition 12.11 such that 
lim I)woo sizc(X>( m) ) = and such that regul(X> (m )) > 0, for all m G N. 
Let p ("0) G H v(m} (Q) d x H vim) (Q) be given by (HU for all m G N. Let us prove 

the existence of a subsequence of (V^ m ') me ^ such that the corresponding sequence 
(^ m ')mEN converges in (L 2 (fl)) 2 to u and the sequence (p^) m m weakly converges 
in (L 2 (fl)) 2 to p, as m — > co. Then the proof is complete thanks to the uniqueness of 
(u,p). 

Using (|3.10j) . we obtain (see [18], [15]) an estimate on the translates of u'™': for 
all m G N, there exists C |g c g^ p , only depending on f2, v, f and g such that 

/H(.T|()- M W)(j))2di < q<n|£|(|£| +4size(X>( m ))), 



mega 



(3.18) 



for k = 1, 



,d, V^G 



transx 



where u( m ' fc ) denotes the fc-th component of v[ m >. We may then apply Kolmogorov's 
theorem, and obtain the existence of a subsequence of (I? < -" i - ) )meN and of u G Hq(£1) 2 
such that (u' m ') m £N converges to u in L 2 (f2) 2 . Thanks to proposition 13.41 we extract 
from this subsequence another one (still denoted u 1 *" 1 * 1 ) such that (p' m ') m gH weakly 
converges to some function p in L 2 (Q). In order to conclude the proof of the conver- 
gence of the scheme, there only remains to prove that (u,p) is the solution of (|3.5|) . 
thanks to the uniqueness of this solution. 

Let tp G (C^(n)) d . Let m e I such that V (m ^ belongs to the above extracted 
subsequence and let (u^- m \p^) be the solution to (|3.7j) with V = . We suppose 
that m is large enough and thus size(2?( m ' 1 ) is small enough to ensure for all K G M 
such that Kn support (y>) ^ 0, then dK n 951 = holds. Let us take v = P V ( m )(p in 
(|3.7[) . Applying proposition 12. 4[ we get 



lim [«<*»> ,P„ (m )d 



Vu(i) : Vy?(x)dx. 
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Moreover, it is clear that 



lim / f(x)-P T>{m) (p(x)dx = I f{x) ■ (p(x)dx, 
n ^°°Jn Jn 



and 



lim 77 / vy n \x) ■ P T ,( m )(p(x)dx = r/ I u(x) ■ ip(x)dx. 
™^°° Jn Jn 

Thanks to the weak convergence of the sequence of approximate pressures, to (|3.1ip 
and to the hypothesis a < 2, we now apply proposition 12. 3i which gives 

(3.19) lim / p {m) {x)div V ( m )(P V ( m) tp)(x)dx = / p(a;)div^(a;)dx. 



The last step is to prove that div(u) = a.e. in f2. Let <p £ C^°(f2) and let m <E 

in d3Zf|). We get p'f^ = 

div vim) (x)(u^)P v{m) ^(x)dx. 



be given. Let us take q = P V ( m )tp in (|3.7|) . We get = — T5 , where 



rp<m) _ 

11 

and 

3^ m) = A size(P^ ) a (pW, P vim w)v- 
On the one hand, the third item of proposition produces 



d „ 

lim T[4 m) = V / ip(x)d l u {l Mx. 
n ^°° m ~1 Jn 



On the other hand, using the Cauchy-Schwarz inequality, we get: 

T| n) < Xsize(V^) a \p^\ v \P v(m) ip\ v 
Therefore, thanks to (|3.11|) and to the regularity of (p (that implies that |Pp(m)y>|p 

rj-'i(m 



remains bounded independently on size(X>(™))) we obtain limn^oo Ti4i"' = 0. This in 



turn implies that: 



i=l 



(3.20) V / if(x)d i u ( - i) (x)dx = 0, for all <p £ C* c °°(n), 



which proves that u £ E(SY). □ 

Remark 3.1 (Strong convergence of the pressure). Note that the proof of the 
strong convergence of p to p is a straightforward consequence of the error estimate 
stated in Proposition \3.6\ below, which holds under additional regularity hypotheses. 



convp 



vandiv 



3.4. An error estimate. We then have the following result, which states an 

error estimate for the scheme (13.71) . 

' * 

Proposition 3.6 (Error estimate in the linear case). Under hypotheses i3.2\) - 
{3.$ . we assume that the weak solution (u,p) of the Stokes problem iS.l]) in the 
sense of definition \3.1\) is such that (u,p) £ H 2 (Q) d x if x (f2). Let X £ (0, +00) and 

13 



a G (0, 2) be given, let T> be an admissible discretization of£l in the sense of definition 
GO and let0>0 such that regul(X>( m )) > 9. Let (u,p) € H v (fl) d x H v (Cl) be the 
solution to \3. 7[ ). Then there exists Cj , which only depends on d, Q, v , rj and 8 such 
that 



(3.21) 



\u-u\\ L2{n) < 



C^e(\, size(2?),p, ', 



eqesterl 



(3.22) 



A sizc(2?) Q \p\l < CM\e(\,size(V),p 



eqester2 



(3.23) 

where 

(3.24) 



\\p-p\\h(Q) < C[^e(A,size(X>),p,u). 



e(\,size(D),p,u) = min (Asizc^)", isizc(P) 2 - Q ) 



x IIP 



H!(n) 



Lff 2 (0) 



eqester3 



epsilon 



Proof. We define (u,p) G Hx>{Q) d x H-p(Q) by u = PpM, which means ma = u(xk) 
for all A' G M. and pk = 777- j K p(x)dx for all A' 6 M. Integrating the first equation 
of (|3~T|) on A G Al gives 



/(x)dx. 



We introduce, for K E M, e 1 ^ = uk — 7777 J K u(x)dx, and, for L 6 Aa: 

R *> L = a^r( fii - - 5^7 Xt V "( a; ) : n A' l( rd7(.T), 

and for a G £ A ' fl £ cxt , #A> = 7^7(0 - "A") - 7777 J a Vtt(x) : iia,^^); 

moreover, we define for L 6 A/a-: e^-^ = ^(pa + Pl) — m ^ JA|i,P( a: ) c ^'( a ')> an d f° r 

(j G £a n fexti £ ct = PA' — 777- J [T p( a; )d7(-i;). Using these notations and the relation 

Effefe nvn^ = 0, we get from (|3.25|) 



77 mj-tiK — f 



E ^AL (PL ~PA\ 

LeM K 



f(x)dx + R K , 



A" 



ith 



R K =i] m K e\- - v I E m A|L-RA,L + E m CT -RA> + £ m <x £ £ n A>- 

We then set du = u — u and dp = p — p. We then get, substracting the first relation of 
the scheme (|3.8p to the above equation, 



(3.26) 



r\ j n Su(x)v(x)dx + v[Su, v]d — J 8p(x)drvT>(v)(x)da 
J n R(x)vdx, Vw G H v {n) d , 



stocontK 



eqester4 
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and, setting v — 5a in (|3.26[) . 



er4 



T) / 6u(x) dx + i/||&t||p 
Jn 



o(x)div-p((5u)(x)da; = / R(x)du(x)dx. 
n Jn 



We now integrate the second equation of (|3.1|) on f e M. This gives 
^ / u(x) • riK iCr d7(x) = 0, V/v G X. 



Using u G i/o (fi) ■ we then obtain 



E A KL ■ (UK + U L )= m K\L£K\L, ^ K G 



with 



2 + 



— / u(x)dry{x) ] • njfi, V/v G 7W, VL G N K - 

K\L Jk\L ) 



We then give, substracting the second relation of the scheme (|3 . 8[) to the above equa- 
tion, 



/ divT>{8u)(x)5p(x)dx = A size(2?) Q (p,p - p)u + T 6 , 
Jo 



with 



Gathering the above results, we get 
(3.27) 



rj J n 6u{x) 2 dx + v\\Su\\l + A size(X>)" \p\l 



A size(X>) Q (p,p)v + Jq R(x) ■ Su(x)dx + 



arerml 



ter5 



Let us study the terms at the right hand side of the above equation. We have, using 
the Young inequality, 



(3.28) 



{p,P)v < -Ml + \P\l < \\P\1 + Cs \\P\\m { n 



eqester6 



We then study J Q R(x) ■ Su(x)dx = T' 7 + T 8 + T 9 , with 



Tfo = 7) J e u {x)-5a{x)dx, 



E ( E ™k\ L Rk,l+ E ^ 



and 



E E m<T £ S n *> ' (5 " A '- 



KEM ae£h 
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Thanks to interpolation results proven in |15| and to (|2.4p . we obtain 



(3.29) 

(3.30) 

and 
(3.31) 



a <C 9 si Z e(2?)i«||^ (n) 
Cio size(X>) 2 ||u|| 2 „ a 



tINI 



2 



J?=(fJ) ■ 4 

isize(2?) 2 ||p||2 rl(n) + 



We then study 2j§) We have = T w - T n with 



which verifies 
(3.32) 

and 



J M = m K\ L £ U K lL (PK -Pl), 

K\LeS iM 



^k\ L s u KI l(pk-Pl), 

K\LeS iM 



2[1T]< U sizc(V) a \p\l + C 13 Uize(Vf- a \\u\\ 2 H2m . 



which verifies 
(3.33) 

;s ji terlObis 

Gathering equations (|3.27p - (|3.33p gives 

Hdujl 2 , + A size(P) Q \ P \ 2 V < qL4${<kimi (v),P 



3r2 . 



where e(A, size(2?),p, u) is defined by ([3~24]) . This in turn yields (l3~2T]l and ((3~22| . We 
then again follow the method used in the proof of Proposition l3.4l Using J n p(x)dx = 
and therefore J Q Sp(x)dx = 0, let v £ i?Q(f2) d be given such that divv(x) = 5p(x) for 
a.e. x £ n and 



ri 2 (n)- 



(3.34) 
We again set 

v® = — [ v {l) {x)dj(x) 7 Wa££, Wi = l,...,d. 
m CT J a 

and we define v £ Hx>(fi) d by 



(») 
,j k 



1 



w (i) (a;)dx, VA' £ M, Mi = 1, 



The same method gives 



11*11 



L 2 (n) 



< 



dp{x)d\YT>{v){x)dx + G^size(V)\p\T>\\v\\ H i 



< / 8p{x)divv{v){x)dx + C 15 sim{Vy\p\ z v + -\\Sp\\l 2 
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eqester7 



eqester8 



eqester9 



eqesterlO 



eqesterlObis 



eqesterll 



We now use v as test function in (|3.26[) . We get 



secf vsnlin 



stnl 



weaksolss 



n 



(x)divv{v)(x)dx = ij I 5u(x)v(x)dx + i/[&t, v]x> + / R(x)vdx. 

Jn Jn 
^^^.^ _____ 

Gathering the two above inequalities, (|3.29[) . (|3.30|) . (|3 . 31 [) and (|3.34|) produces 



1*11^(0) < dlfclli 



Applying 



(fi) ■ 
+C 17 ||&|||,+ 

and (|3~3-| gives (fBT^j) . □ 



Ci 6 size(D) [\\p\\ H i m 



\H 2 (n) 



*(V?\ P \l. 



Remark 3.2. In the above result, it suffices to let a = 1 to obtain the proof of an 
order 1/2 for the convergence of the scheme. We recall that this result is not sharp, 
and that the numerical results show a much better order of convergence. 

4. The finite volume scheme for the Navier- Stokes equations. Before 
handling the transient nonlinear case, we first address in the following section the 
steady-state case. 

4.1. The steady-state case. For the following continuous equations, 

d 

;) + dip + ^TuWdjii® = /« in n, for i = 1, . . . , d, 

3=1 



r^vy' — v 



(4.1) 



divti = 2_,®i u = m ^- 



with a homogeneous Dirichlet boundary condition, we define the following weak sense. 

Definition 4.1 (Weak solution for the steady Navier-Stokes equations). Under 
hypotheses H3.ty )- pT4\ ), let E(fl) be defined by U.8\) . Then (u,p) is called a weak 
solution of 0771) if 



(4.2)1 



u £ E(Q), p G L 2 (n) with [ p(x)dx = 0, 

Jn 

•q I u(x) ■ v(x)dx + V I Wu(x) : Vv(x)dx 
Jn 

p(x)divv(x)dx + b(u,u,v) = / f(x) ■ v(x)dx Vw G Hl(Vt) d , 



where the trilinear form b(., ., .) is defined by \1.10f) . 

We now give the finite volume scheme for this problem. Under hypotheses (|3.2 
(13.41) , let T> be an admissible discretization of SI in the sense of Definition 12.11 We 



nstocontss 



nstocontf ss 



mega 



introduce Bernoulli's pressure p + \u 2 instead of p, again denoted by p, and for any 
real value A > and a G (0, 2), we look for (u,p) such that 

(u,p) G H v {ri) d x H V (Q) with / p{x)dx = 0, 



(4.3) ^ 



?7 / u(x) ■ v(x)dx + v\u, v]x> + — I u(x) dbfT>{v)(x)dx 



p(x)d\VT>{v){x)dx + br>(u, u, v) 



d\vxi{u)(x)q(x)dx = — A size(2?) Q {p,q)r> 
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f{x)-v{x)dx VveH v {ny 
Vq G H v (n) 



schvfnlss 



where, for u,v,w G Hx>{£l), we define the following approximation for b(u, v, w) 

(4.4) br>{u,v,w) = ^ ^ ^ {A KL ■ (u K + u L j) ((v L - v K ) ■ w K ) 

rem LeM K 

iRr.h vf nl ss 

System (|4.3|) is equivalent to finding the family of vectors {uk)kgm C K , and 
scalars (pk)kem C K solution of the system of equations obtained by writing for 
each control volume K of M.: 



-j—(u l -uk)-v ^ 



(4.5) { 



^ (Arx • (-(% + ml))) {ul - uk) 

L^Mk 

A kl (Pl-Pk) - ~ A KL ( 



Iff,! 



-(0-ttx) 



LeM K 



LeM K 



) = / f(x)da 



K 



A KL -{u K + u L )-Xsizc{V) a 



LeM K 

supplemented by the relation: 



LeM K 



™-K\L 
&K\L 



(j>l -Pk) = 



mjf pa- = 



Defining p K = Pk - u 2 K /2 and p CT = (p K + Pl)/2 if a = K\L, p a = p K if 
a G £ ext n and using the fact that J2cre£ K m c n K,<7 = 0, one again notices that: 
J2l£AT k Akl (pl ~Pk) is in fact equal to J^ <Te£l( m (T p (T nA>, thus yielding a conserva- 
tive form for the fifth and sixth terms of the left handside of the discrete momentum 
equation in (|4.5|) . Defining u a = [uk + ul)/2 if er = K\L, u a = if a G £ e xtn£ir, one 
obtains that the nonlinear convective term YIlgNk (Akl ■ {\{ u k {ul~ uk) is 

equal to J2ae£ m cr( n K,tr ■ u a )u a — mKUK(drvx>u)K', one may note that (divpu)^ = 
^2 a £ £ m CT nA',cr • u a . Hence the nonlinear convective term is the sum of a conserva- 
tive form and a source term due to the stabilization (this source term vanishes for a 
discrete divergence free function u). 

Let us then study some properties of the trilincar form 6-p. First note that the 
quantity bx>{u,v,w) also writes 

(4.6) bv(u,v,w) = - Y {Akl ■ {uk + u l )) ((v l - v k ) ■ (w l +w k )) 

^ K\L££ int 

We thus get that, for all u,v G Hj)(D,) d , 



br>{u, v, v) 



\ E ^KL 



{u K + U L )){{v L f -{v K ?) 



(4.7) 



K\Le£ h 

1 

2 



v{x) 2 divx>(u)(x) da; 



We get in particular, that, for all u G Ev(Cl), bx>{u,u,u) — 0, which is the discrete 
equivalent of the continuous property 
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def tridc 



schvfNSS 



def tridebis 



superb 



Remark 4.1. [Upstream weighting versions of the scheme] All the results 
of this paper are available, setting Fkl( u ) = ^kl ' { U K + u l) an d considering, for 
u,v,w G Hu(n), 

b^ s (u,v,w) =br>(u,v,w) + - ®kl\Fkl{u)\ {v l - v K ) ■ (w L - w K ), 

K\Le£ int 

with, for example, &kl = max(l — 2v ™'^ /\Fkl{u)\, 0). We then get, for all u, v € 
Ht>(Q), the ineguality 

6p PS (w, v, v) > — — I v(x) 2 divi>(u)(x)dx, 
2 Jn 

which is sufficient to get all the estimates of this paper, together with the convergence 
properties of the scheme. The use of such a local upwinding technique may be useful to 
avoid the development of nonphysical oscillations only where meshes are too coarse. 

The following technical estimates are crucial to prove the convergence properties 
of the scheme. 

Lemma 4.2 (Estimates on &£>(., ., .) by discrete Sobolcv norms). Under hypotheses 
U ,3\ )- [L7\ ), letT> be an admissible discretization of fix (0,T) in the sense of definition 
\4-8[ and 6 > such that regul(D) > 8. Then there exists ( .x >!() and C'm > 0, only 
depending on d, 9 and £1, such that 

(4.8) bv(u,v,w) < qj^ \\u\\ L 4 {n)d \\v\\ v \\w\\ L 4 {n)d < qjT?] \\u\\ v \\v\\v 1Mb- 



Proof. The quantity bx>{u,v,w) reads 



bv(u, v,w) = j E ( Wk ' ( Wl _ VK ^ ^ A |L ~ x k) ■ ( u k + u L )) 



4 ^ , d K \ T 



Applying the Cauchy-Schwarz inequality twice and using the fact that {xl — xk) 2 

;LeA/ K d KIL KL 



d 2 KL and that, for any admissible discretization Xl<ea/ TT^ d 2 KL < d yield: 



b v {u,v,wf < C 2Q V V -^( WK f( XL - XK )\2{u K ) 2 + 2{u L ) 2 ) 



( E E t^(-) 2 ( 

(E E ^ L -vkA 

< c 2 i I ^2 m K\w K \ i j [ 



1/2 / \ 1/2 

4 I 



m K \u K \ 
kem J 

The inequality (|4.8[) is now a straightforward consequence of the following discrete 
Sobolcv inequality, which holds under the same regularity assumptions on the mesh 
(see proof in [10] or [El pp. 790-791]): 

(4.9) IMUo(n) < C22 1Mb- 
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degtop 



estlll2ss 



Remark 4.2 (Two dimensional case). In the case d — 2, it may be proven setting 
■■ 2,p = p' = 2 in the proof pi 91 of £f5|/ ; that 



and therefore, that there exists C24 > j), only depending on d and 0, such that 
bT>{u,v,w) < C^24\\\v\\v (\\u\\y ||u|| L 2 (n) \\w\\t> \\w\\ L 2 {n) ) 1/2 . 



1 1/2 



This is a discrete analogue to the classical continuous estimate on the trilinear form. 

The existence of a solution to the scheme (|4.3|) is obtained through a so-called 
"topological degree" argument. For the sake of completeness, we recall this argument 
(which was first used for numerical schemes in [T7]) in the finite dimensional case in 
the following theorem and refer to [TT| for the general case. 

Theorem 4.3 (Application of the topological degree, finite dimensional case). 
Let V be a finite dimensional vector space on R and g be a continuous function from 
V to V . Let us assume that there exists a continuous function F from V x [0, 1] to V 
satisfying: 

1. F{-, 1) = g, F(-,0) is an affine function. 

2. There exists R > 0, such that for any (v,p) G V x [0, 1], if F(v,p) = 0, then 
\\v\\ v £R. 

3. The equation F(v,0) = has a solution v £ V such that \\v\\v < R- 
Then there exists at least a solution v £ V such that g(v) = and \\v\\v < R- 
Here g(v) = represents the nonlinear system (|4.3p . and we are now going to 

construct the function F and show the required estimates. Note that here, the use of 
Bcrnouilli's pressure leads to simpler calculations. 

Proposition 4.4 (Discrete Hq (fl) estimate on the velocities). Under hypotheses 
i3.2\ )- [3~4\ ), letT) be an admissible discretization of fix (0,T) in the sense of definition 
\4-S\ Let A G (0, +00) and a £ (0, 2) be given. Let p £ [0, 1] be given and let (u,p) € 
(-ffp(Sl)) d x Hj)(fi), be a solution to the following system of equations (which reduces 

to jr ' 



as p = 1 and to \3. H\j as p = 0) 
(u,p) G H v (n) d x H V (Q) with I p{x)dx = 0, 



77 / u(x) ■ v(x)dx + v[u, v]d + — / u(x) 2 divxi(v)(x)dx 



+p bT>{u,u,v) — / p(x)divT>(v)(x)dx = / f(x) ■ v(x)dx 



divx>(u)(x)q(x)dx = —A sizc(P)" (p,q)i 



Vi> g H v (n) d 

\lq G ffu(n) 



(4.10) 

Then u and p satisfy the following estimates, which are the same inequalities as ob- 
tained in the linear case (inequalities \3.10\l and &3.11\) ): 

v\\u\\t> < diam(n)||/|| (L 2 (n)) d 



schvfnlssrho 



v A size(V) a \p\l < diam(0) 2 ||/|| 2 L2 
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Proof. The proof is similar to that of Proposition [321 using the property (|4.7p on 
the discrete trilincar form. □ 



2 

f{ x )-v(x)dx \/veH v (ny 



We are now in position to prove the existence of at least one solution to scheme 

631). 

fimega 



Proposition 4.5 (Existence of a discrete solution). Under hypotheses US. 2\) 
j3.4\) . let D be an admissible discretization of x (0, T) in the sense of definition 
\4-8\ Let A G (0, +oo) and a G (0.2) be given. Then there exists at least one (u,p) £ 
(H v (ty) d x H v {fl), solution to ( fpf ), 

Proo/. Let us define V = {(u,p) G (H v (Q)) d x i?p(0) s.t. J n p(x)dx = 0}. 
Consider the continuous application F : V x [0, 1] — > V such that, for a given 
(u,p) € and p G [0, 1], (u,p) = F(u,p,p) is defined by 

u(x) ■ v(x)dx = T] I u(x) ■ v(x)dx + v[u,v]x> — / p(x)div£>(t;)(x)dx 
Jsi Jo 

+P ( - /" u(x) 2 divp(v)(x)dx + 6x>(u, u, u) 
n 



p(x) • q(x)dx = / divx>(u)(x)g(x)dx + A size(2?) a (p, g)p Vg G ifc>(f2). 

It is easily checked that the two above relations define a one to one function F(., ., .). 
Indeed, the value of u K ^ and px for a given K G .M and i = 1, . . . , d are readily 
obtained by setting «W = . uU) = for j ^ i, and g = Ik- 

The application i r (., ., .) is continuous, and, for a given (u,p) such that F{u,p, p) = 
(0,0), we can apply proposition 14.41 and (|2.5j) . which prove that (u,p) is bounded 
independently on p. Since F(u,p,0) is an affine function of (u,p) (indeed invertiblc, 
see corollary [231) 7 we may apply Theorem l4.3l and conclude to the existence of at least 
one solution (u,p) to (|4.3[) . U 

We then have the following strong estimate on the pressures. 

PROPOSITION 4.6 (L 2 estimate on pressures). Under hypotheses \3.$ty - \34\ ), let 
T> be an admissible discretization of in the sense of definition \2.1l and let 9 > 
such that regul(P) > 6. Let A G (0, +oo) and a G (0,2) be given. Let (u,p) G 
HT>(fl) d x Hx>(U) be a solution to J^.3[ ). Then there exists C25 , only depending on d, 
f2, r\, v , A, a and 9, and not on size(2?) ; such that the following inequality holds: 



(4- 11 ) IIp!Il 2 (J2) < (ll/|[(L 2 (fi)) d + (ll/ll(L2(J2))<i) 2 ) 



Proof. We may follow the proof of proposition 13.41 until (|3.17[) , which is changed 



to: 



(4.12) 



p(x)div£>(z;)(x)dx = r\ I u(x) ■ v(x)dx + v[u, v\d — I f{x)-v(x)dx 



+ 2 / W (' T ) divv(v)(x)dx + bv(u,u,v). 



vdansschvfnlss 
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We again apply the discrete Poincare inequality (|2.4[) . (|3.15p . (|3 . and we use 
We get the existence of CEg^toij^ipbending on d, O, /, rj, v, A and 6, such that 

Iblli^n) - q^size(V)\p\ v C^\\p\\ L 2 {n) < 

qsi (Nb + ll/IUw + 1Mb) 

We now apply (|3.10p and (|3.11[) . which yields the conclusion. □ 

lyf nlss 



We now can state the convergence of Scheme (|4,3|) . 

Theorem 4.7 (Convergence of the scheme). Under hypotheses i3.ty) - pT^\ ), let 
(•p( m )) mgN be a sequence of admissible discretizations of Q in the sense of definition 
\2.1l such that size (r>( m )) tends to as m — > oo and such that there exists 9 > 
with regul(2? (m) ) > 6, for all m G N. £e£ A G (0, +oo) and a G (0,2) be given. 
Let, for all m G N, (u( m \p( m >) G {H v(m) (n)) d x H v{m) {VL), be a solution to 



with T> = T>( ' \ ls Then there exists a weak solution (u,p) of \4-.l\ ) in the sense of 
definition \4-l\ and a subsequence of (f^)meN) again denoted (£>( m )) meN; such that 
the corresponding subsequence of solutions (u^) m gN converges to u in L 2 (£l) and 
(pi™) — i(u( m )) 2 )„ lSN weakly converges to p in L 2 (f2). 

Proof. Since the same estimates as in the linear case are available in the steady 
nonlinear case, the proof of proposition 13.51 holds for all the terms of (|4.2|) which are 
present in (|3.5j) . We only have to prove that for a given <p G (C£°(f2)) d , as m — > +oo: 

T{™' = \ u^ m \x) 2 &Yv V { m ) (P T ,( m )(p)(x)dx tends to / u(x) 2 divip(x)dx 

Jq Jn 

and 

= b v (ut m \ U ( m \P v(m) tp) tends to b(u, u, tp). 



Thanks to the convergence in L 2 (n) of (u (m) ) me 



to u and to the discrete Sobolev 
v \\li(Q) < C 2 7<jMh^ for au v e H V ( m) {Vt) and all g < 6 (see [HI p. 
using (|3.10p the convergence in L 2 (fl) of ((M^ m ^) 2 ) me N to u 2 . We now 
1, . . . ,d, the sequence (Ppim) tp® ) m eN satisfies the hypotheses of 
Hence, Vd(»o P©cm) weakly converges to V^' 1 * 1 in L 2 (Q) d . One 
has divx.w = EiLi v£V*> for all u G (H v (n)) d such that it K = if £jr n £ cxt ^ 0. 
Hence divp(m) (Pptm) weakly converges to div<^ in L 2 (r2), thus providing the limit 



inequalities 
790]), we g( 
remark that for i 
Proposition [5 



of Tr 



(m) 



d 

J dg.fjtridcbis , , 

Thanks to (|4.6[) . setting for simplicity I? = T^™- 1 , we have: 



bv(u,u,P v ip) = I f^s^) - Ij^TJ 



with: 



Tr 



Ml 



EE 

fc=li=l 



, ■ u k ){{ul - «jf) • <p(xk)) 
(x)V$ (w (fc) ) {x)P v ^ k) (x)dx 



E 



(Ar- l ■ (u L - u K ))((u L - u K ) ■ (<p(xk) ~ V^i))) 
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Thanks to the convergence in L 2 {Cl) of (w^P-po) <^) m(E N to uip, we get from propo- 
sition E5D that: 



d c 
fc=li=l J " 



u w (x)$uW (x)<^ w (z)da; = &(u, u, if). 



We have: 



/ 



[El 



E 



m A '|L 



and therefore, since \(p(xk) — V>(xl)\ < dKLC v sizc(T>) where C v is a bound of V<p in 
L°°(Cl) d , and since cZrx < 2size(X>), the following estimate holds: 



(, 

m 



< 4size(P) 2 C ¥ 



\T>- 



Therefore, (j3~TUj) yields: 



lim Tr 



(m) 



El 



o. 



which concludes the proof of convergence. □ 

4.2. The transient case. We now turn to the study of the finite volume scheme 
for the transient Navier-Stokes equations, the weak formulation of which is given in 

(HH). 

We first give the definition of an admissible discretization for a space-time domain. 

Definition 4.8 (Admissible discretization, transient case). Let CI be an open 
bounded polygonal (polyhedral if d = 3) subset of M. d , and dCl = Cl \ Cl its boundary, 
and let T > 0. An admissible finite volume discretization of CI x (0, T), denoted by 
T> , is given by T> = (A4, £, V, N), where (A4,£,P) is an admissible discretization of 
CI in the sense of definition \2.1\ and jVgNj is given. We then define St = T/N , and 
we denote by size(2?) = max(size(A4, £, V), St) and regul(P) = rcgul(A^, £, V). 

Under hypotheses (|1.3D - (11.7p , let T> be an admissible discretization of CI x (0,T) 
in the sense of definition 14.81 and let A E (0, +oo) and a <G (0, 2) be given. We write a 
Crank-Nicholson scheme for the time discretization, and follow the nonlinear steady- 
state case for the space discretization; the finite volume scheme for the approximation 
of the solution (fTTTI) (fTT2]) is then: 



(4.13) 



u Q G H v (Cl) d , 

uo k = / u; n i(x)dx, V-FT G M, 

ms- Jk 



and, again using Bernoulli's pressure p + \u 2 instead of p, again denoted by p, 
(u n+1 ,p n+ i) G (H v (Cl)) d x Ht>(Q), 

J n p n+ l{x)dx = 0, U n+ L = \(u n+ x +Un), 

(u n+ i(x) - u n (x)) ■ v{x)dx + vSt[u n+ i ,v]t> 

St 

(4.14) 



-St / p n+ i (x)diw(v)(x)dx + - / u n+ i(x) divD(v)(x)dx 



+Stb v (u n+ i,u n+ i,v) 



2 

(n+l)St 



f(x, t) ■ v(x)dxdt. 



divT>(u n+ i)(x)q(x)dx = -A size(X>) a (p n+ x,q)v, 
Vu G H v (Cl) d ,Vq G H v (Cl), Vn G N. 
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In (|4.14p . wc consider the approximation of 6-p given by (|4.4[) . We then define the set 
Hx>{Q x (0, T)) of piecewise constant functions in each K x (ndt, (n + l)St), K G A4, 
n G N, and we define (u,p) G i?x»(f2 x (0, T)) by 

u(x, i) = u n+ 1 (x), and p(x,t) = p n +^(x), for a.e. (x, t) G fix (n$, Vn € N.| 

(4.15) 

Remark 4.3 (Time discretization). It is wellknown that the Crank-Nicholson 
discretization is implicit. If we use the 8 scheme: u n +± = ^ u n+i + (1 — 0)u n , with 
G [3,1], the convergence proof which follows applies with a few minor changes. 
Variable time steps may also be considered. 

Let us now prove the existence of at least one solution to scheme (|4. 13() - f|4. 15|) . 
Proposition 4.9 (Existence of a discrete solution). Under hypotheses U.3\) - 



. let T> be an admissible discretization of fl x (0, T) in the sense of Definition 
iLet A G (0, +00) and a G (0,2) be given. Then there exists at least one (u,p) € 
{Hv(n x (0,T))) d x H V (Q x (0,T)), solution to £p?l\ )- [p]fy . 

Proof. We remark that, for a given n = 0, . . . , N — 1, taking as unknown u n+ i, 
and noting that u n+ \ = 2u n+ i — u n , Scheme (|4.14j) is under the same form as scheme 
(|4.3p . with r\ — \ and with a term in u n included in the right hand side. Therefore 
the existence of at least one solution follows from proposition 14.51 □ 
We then have the following estimate. 

Proposition 4.10 (Discrete L 2 (0,T; Hq(Q,)) estimate on velocities). Under hy- 
potheses A1.3\) - p77\ ), let T> be an admissible discretization of SI x (0,T) in the sense 
of definitional^ Let A £ (0, +00) and a £ (0,2). Let (u,p) G (H V (Q x (0,T))) d x 
Ht>(£1 x (0,T)), 6e a solution to J^.j3| )- f^.i5| ). TTien t/iere exists C29, > 0, onZy 
depending on d, f2 ; j/, Uq, /, T suc/i i/ia£ the following inequalities hold 



(4.16) 

(4.17) 

and 

(4.18) 



MI,L^(UT;L 2 (Q) d ) < 
MU 2 (0,T;ff 13 (n)< 1 ) < 



N-l „ T 

A size(X>) tt £ Sb\ Pn+ i \% = A size(P) Q / |p(-, df < 



estimutinf v 



estimut 



estimpt 



Proof. Let p = 1, . . . , N. We get, setting v = w n+ i in the first equation of (|4.14|) . 
summing on K G M. and n = 0, ... ,p — 1 in the first equation of (|4.14[) and using 
property (|4.T|) . 

p— 1 ^ p— 1 

^„(x) 2 )dx + v (S[l 



1 p_1 /• 

p-1 



n+1 



(x) 5 



n=0 
p-1 



V« / p n+ i(x)divx>(u n+ i)(x)<£x = V/ / /(x,i) • u„ + i(x)dxdi, 

This leads, setting g = p n+ i in the second equation of (|4.14[) . to 

(u p (x) 2 



(4.19) 



p-i 

Mo(x) 2 )dx + V &[ u n+l 
P- 1 r p<S 



,u n+ i jx>- 



\size(vrY,*\Pn+±\ 2 v= \ 



/(x, t) ■ it(x, t)dxdt. 



meqen 
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Setting p = N in (|4.19j) gives (|4.17[) and (|4.18[) . The discrete Poincare inequality 
and the inequality ||ito|U 2 (n) d < \\ u im\\L 2 (n) d g iv e 



dis 



L P\\L 2 (Sl) d 



< diam(O) 2 



2;/ 



L 2 (fix(0,T)) d 



i n+ i \\L 2 (ay 



< 



\Uim\\ L 2 (Q) d, Vp=l,...,N, 
\ L 2 {n)d + H^n+ilU^n)") for all n 



which proves (|4.16|) . since 
0,...,N-1. Q 

We then have the following estimates on translations. 

Proposition 4.11 (Space and time translate estimates). Under hypotheses 
il.3\ )- [T7?\ ), let T> be an admissible discretization of Q x (0,T) in the sense of def- 
inition Let X G (0, +oo), a G (0,2) and 8 > 0, such that regulfD) > 9. Let 



(u,p) G [H v {n x (0,T))) d x H v (fl x (0,T)), be a solution to [Jl3\ )- [4l5^ . We de- 
note by u the prolongment in M. d xl of u by outside of fix (0, T). Then there exists 
C29 > and C30 > 0, only depending on d, £1, v, A, a, uq 7 f , 6 and T such that the 
following inequalities hold: 



(4.20) 

and 

(4.21) 



\H- + t-)-u\\ 



< 



|£|(|£|+4size(M)), V£g 



\\u(-,-+t) - u\\ L i { 



-L 1 



< 



1 1/2 



VrG 



estimtrsp 



estimtrti 



Proof. In the following proof, we denote by C{, where i is an integer, various 
positive real numbers which can only depend on d, f2, is, A, a, uq, /, and T. 
Inequality ([420]) is obtained from (|47TT]) (see [l5j). Let us prove (|4~2Tj) . Let r € (0, T) 
be given. We define the following norms on (i?p(f2)) rf , by: 



Vwg {H v (n)) d , 



(4.22) 

and 
(4.23) 



23, A 



Asize(X>) Q 



sup 



divT,(w)(x)q(x)dx, q G i?x»(fi), \q 



v 



V w G (£Mft)) d , 

||tu||*,z>,A = sup { J n w(x) ■ v(x)dx, v G (HT>{Q)) d , |M|x>,a = l} • 



We then have, for a.e. t G (0,T), 

\\u(-,t + r) - u(-,t)\\ 2 L2{n)d < \\u(-,t + r) -u(-,i)||i>,A||u(-,i + r) - u(-, t)||*,x>,A, 
and therefore, thanks to the Young formula, 

( ] +^ 7 \\u(.,t + T)-u(;t)\U,V,X- 



teps 



We get, from (jL14| . for all q G i?x>(fi) and for a.e. t G (0,T), 



normeel 



divT>(u(- ,t))(x)q(x)dx = -A sizc(X>) Q (p(-,i), g)x>, 

25 



which proves, using (|4.22[) . that 

R>t)llz>,A< \H;t)f v + Xshe(V) a \ P (;t)\ 2 v . 
Using the Cauchy-Schwarz inequality, we have that: 

t-t \ 2 pT 

|K-,t + T-) -u(.,t)||i,,xdtj <4Tjf \W>t)\\l,\dt, 

, ,.mt last, i nipt. 

and therefore, using (|4.17|l . (|4.18|) . 



(4.25) 



T-T 



|K,t + T)-u(-,t)||i ? ,Adt<C , 31. 



We now study + r) — i)||±,z>,A- We can write, for a.e. is (0,T — r) and 

x g o, 



JV-l 



(a;,i + r) -w(x,<) = i E T ) + Xn+i(*ji"))(«n+i( a; ) -«n(^)), 



bordlambd 



where, for all n € N and t G (0,T), x„(i, r) = 1 if n& € [i,t + r[, and x„(t, r) = 
otherwise. This implies 



(4.26) 



IK-,t + r)-u(-,t)||*,i>,A < 



hT,n=o(Xn(t,T) +Xn+l(t,T))\\u n+ i - U„||*,D,A- 



Let us then obtain a bound for ||u n +i — ttn||*,z>,,V Using the scheme (|4.14[) . we get 
that, for all v G (H v {fl)) d , 



(u n+ i{x) - m„(j;)) • u(a;)dx 



(n+l)& 



7 / (V/ 



/(x, t) • v{x)dxdt 



—v8b{u n+ i ,v]t> + & / p n+ i(x)divx>(v)(x)dx 
2 Jn 2 

St f 

-7 / u 2 i+L div v (v)(x)dx - Stb v (u n+ i,u n+ i,v). 



(4.27) 

Using the definition of divp, the fact that J2aes K m cr n K,u = 0, and the Cauchy- 
Schwarz inequality, there exists C32 such that: 



u 2 n+ i(x)div v {v)(x)dx < <^\\u 2 n+ i\\ L 2 (n) \\v\\ v . 



The discrete Sobolev inequality (|4.9[) leads to 



ght 



H< + ilU=(n) < ElK*2-*) a H^(n) = Ell«2-*Hi*(n) ^ C » K+l^ 



i=i 



We take |M|x>,a = 1 and note that, from Definition (|4.22j) . we obtain that \\v\\t> < 1, 
and that f n p n+ i(x)dTVx>(v)(x)dx < (A sizc(I?) Q ) 1 ^ 2 |j>„ + i We then pass to the 

26 



supremum in (|4.27|) . Using the Cauchy-Schwarz inequality, the discrete Poincare 
inequality, and (|4.8[) . this yields: 



\\u n +i - u n \\^ v ,x < V&diam(fi)||/|| i 2 (0x(na ( „ +1)lS)) (S 

+Ai/||« B+ 4|| D + (A sizc(P)«) 1/2 \ Pn+1 \ v 

l ... 

+*(2 C l32] <: l33l + C M)W u n+^ He- 



Summing the above equation for n = to N — 1, applying the Cauchy-Schwarz 
inequality to all terms of the right hand side except the last, using (|4.17p and (|4. 18[) . 
we get that there exists C34 such that 

JV-l 



\\u n +i - u n ||*,D,A < Q34l- 



Hence, noting that for all n = 0, . . . , N, L T Xn(t, r)dt < r, we have: 
2 



1 j.T-tN-1 

9 / 51 (X"^' 7 ") +Xn+l(*' r ))ll'"™+l -"«ll*,X>,Adt < Cj5|]T, 



which proves, using (|4.26j) . 

,T-r 

(4.28) / \\u(;t + T)-u(;t)\U, V , X dt< 



Thanks to P~2"4"]l . P~2"5|) and P~2"g]l . we obtain that 



T-r 



||«(-,t + r) - u(-,t)|| i 2( U) ddi < C 3B ftff . 



Using (|4.16p . we have 



\H;t + T)-u(-,t)\\ L 2 (Q)d dt= [ \\-u(-,t)\\ L 2 {ny dt<QmT< vWrCjm 

T-r JT-t 

and a similar inequality holds for f \\u(-,t + r) — fi(-, t)\\ L 2^ddt. This thus gives 
(I4.2ip . for any t 6 (0, T). The case r > T is obtained again using (|4.16|) . and the case 
t < is obtained from t > by the change of variable s = t + t. This completes the 
proof of P~2T|) . 

Theorem 4.12 (Convergence of the scheme). Under hypotheses il.3\) - p77\ ), 
let 9 > be given and let (2?^ m ^)meN be a sequence of admissible discretizations of 
n x (0,T) in the sense of definition \JH[ such that regul(p( m )) > and size(r>( m )) 
tends to as m — > 00. Lei A £ (0, +00) and a <E (0, 2) fee given. Let, for all m € N, 
(«W,j)( m ») e (H v(m) (n x (0,T))) rf x H v(m) (n x (0,T)), fee a so/wiion fo g^f- 
J^._/5[ ) T> = T>( m > . Then there exists a subsequence of {T>^ m ^) m< zm, again denoted 
(2?(" l )) mgNj such that the corresponding subsequence of solutions {iS m ^) m ^ converges 
in L 2 (Q, x (0,T)) to a weak solution u of \ l.l\) - n~Jfy in the sense of definition \l.l\ 

Proof. Let us assume the hypotheses of the theorem. Using translate estimates 
(|4.20p and (|4.21j) in the space L 1 (M. d xl). we can apply Kolmogorov's theorem. We get 
that there exists u € L 1 (f2 x (0,T)) and a subsequence of (X>( m )) me jsj, again denoted 
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(£>( m )) m6Nj s U ch that the corresponding subsequence of solutions (u' m ') ra gN converges 
in L\n x (0,T)) to u as to -> oo. Using BT71). we get ||u (m) ||L2(0,T ; H Om (J2)) < ^28], 
forallme N, which gives, using the discrete Sobolev inequalities, ||u( m '||.i,i(o,T;.L 4 (n)) <| 
Cbebnft'i' all m E N. Using a classical result on spaces L p (0, T; L g (0)), we get that 
(«' m ')mgH converges in L 1 (0,T;L 2 (fi)) to u as ;» — > do. Thanks to (14.161) . we have 
ll" (m) IU~(o,T;L2(n)<*) < for all m G N. The same result on spaces £ p (0, T; 
implies that (u( m )) m(E N converges in L 2 (0, T; L 2 (0)) to u as to — > oo. We can there- 
fore pass to the limit in (|4.20p . The resulting inequality implies u G L 2 (Q, T; Hq (fl) d ) 
(see US]). Passing to the limit in (f4~T6]l leads to u € L°°(0, T\ L 2 {VL) d ). 

Lelj. us now prove that u is a weak solution of l|l.l|l - (|1.2|) in the sense of definition 

o 

Let if G C^°(fi x (-oo,T)) d be given, with divv3(x,i) = for all (x,t) G x 
(—oo, T). Let 2?( m ' be a given admissible discretization extracted from the considered 
subsequence. Omitting some of the indices to for the simplicity of notation, we then 
set v = Pd (/?(•, n5t) in (|4.14|) . and we sum for n = 0, . . . , N — 1. We thus get 

(4.29) T^g' + Sjgj* + T^fg + /[jTT] + ^ = T^jf, 



wi 



th 



iV-1 . 

= ^2 ( u n+i{x) - u n (x)) ■ P V (p(x,nSt)dx, 

n=0 ^ n 



N-l 



T 17~' = H Ak+i'^^''"*)]^^ 



n=0 



T- 



18 



iV-1 . 

^2 * / Pn+i( a; ) div x)(fx)'^(-, n#))(x)da 



1 JV_1 /" 

2j^ = 2 X! * / w n +i(^) 2 divD(PD^(-,n«))(a;)da 

n=0 



AT-1 



2^f = ^ &&x>(u n+ i,iV4' P ^('' n< *))' 



and 



n=0 



ri J ndt J Q 



/ f(x,t) ■ Pr><p(x,nSt)dxdt. 
Jn 



In the following, we denote by Cj various positive reals which can only depend on d, 
f2, T, itjni, /, i/, 9 and A. We first start with the study of T$yj\ We classically have 



(see 



(4.30) 



lilll jry 



o Jn 



Vu(x,i) : V<p(x, t)dxdt. 



cvt2 



28 



The proof that 

Jp^tiijipt. lypakcOTivgrarl 

is a consequence of ()4.18|1 and of a direct adaptation of Proposition 12.31 to time- 
dependent functions. Let us now prove that 

(4.32) lim TiS ) = 0. 

Since (u*-™* 1 ) 2 tend to u 2 as m — > oo in L 1 (fi x (0, T)), the same argument as in the 
steady state case (see proof of theorem 14. 7| provides (|4.32jl . 

We now turn to the study of 2 j2Q} Following the proof of proposition 14. 7\ the 
proof that 

(4.33) I™oc J M = j Ku(;t),u(;t),<p(;t))dt, 

is a direct consequence of the convergence of u to u in L 2 (il x (0, T)) and proposition 
12.31 The study of T ^j\ is classical, and we have 

(4.34) lim r 4m = [ [ f(x,t) ■(p(x,t)dxdt. 
Let us now prove that 

(4.35) lim Trfp\ = — / / u(x, t)dt(p(x, i)dxdi — / Ui n i(x)ip(x,0)dx. 



Indeed, we have 



with 



jv-i . 

" ^ / u n+ i (x) • {Pt>p(x, (n + 1)4) - Pv^{x, nSt))dx. 

n=0 *' n 



and 



JV-l 



= y2 (u n +i(x) - u n (x)) ■ (P v <p(x, (n + l)St)- PD<p(x,n&))da 

We classically have 

lim / ito(x) • Ppf{x, 0)dx = / Wi n i(x)</?(x, 0)dx. 
m - > °° Jo Jo 

We also easily have, thanks to the convergence properties of u^ m \ that 



m — 'oo 



lim = / / u(x,t)dt<p(x,t)dxdt. 

' Jo Jn 
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Let us prove that the term T^j ' 1 tends to as m — > oo. We have = ^i-T' — 
with 



iv-i r 

] = ^2 (u n+ i(x) ~ u n (x)) ■ Pv^{x, (n + l)dt)dx. 

n=0 ^° 



Thanks to the limits given by (|OD|) . (|QTj) , P~53"| and (|Q4"|) . and thanks to 

(|429|) , we obtain that lim TrfS ) = T25, with 

7 &B= - U Y. Wu^{x,t)-W^{x,t)6xdt- b( U (;t),u(;t),tp(;t))dt+ 

■ , Jo Jn Jo 



i=l 
T 



Jn 

f(x) ■ <p(x, t)dxdt. 

Jo Jn 

Irvt? Irvt.3 IrvtSh \r\rtA Irvt.R 

Since (|430]) . (|4~3Tj) . (j4~32]) . (|433]l and (|434|) are available as well, replacing Ppipi^nSt) 
by Px>^(-, (n + 1)&) in TQj], TJgj and Tj^jj, we also get using <|4.14|) with 

v = Pr>(p(-, (n + I) St) that lim = Tf^g) . Thus we get that lim m _,oo = 0, 

which concludes the proof of (|4~55l) . Thanks to (|4T29]) . (j435j) . (|QU| . (|4~31"]l . (|4~32"1) . 
()4.33p and (|4.34|) . we thus obtain (| 1 . 9|) . provided that we can prove 

divu(x, t) = 0, for a.e. (x, t) G ft x (0, T). 

This last relation can be shown, following the proof of (|3.20[) . This completes the 
proof of the above theorem. □ 

Remark 4.4. Using the above proof of convergence, we get the energy inequality 
for d — 2 or 3 from inequality J^. 19\ l, since we have the property 



/ / (W%,t)) 2 dxdi<liminf V &[u%t«%$U m 
Jo Jn " wo ° ^ +5 + ^ 



5. Numerical results. An industrial implementation of a colocated finite vol- 
ume scheme may be found in [1] for instance, where complex applications are consid- 
ered. Focusing in this paper on properties of convergence and error estimates, some 
simple numerical experiments are described here to observe the convergence rate of 
Schemes (|3. T[) and (|4.13|) - (|4.14|l with respect to the space and time discretizations. To 
that purpose, we use a prototype code where the nonlinear equations arc solved by 
an underrelaxed Newton method, and the linear systems by a direct band Gaussian 
elimination solver. This code handles Stokes or Navicr-Stokcs problems with various 
boundary conditions, using non uniform rectangular or triangular meshes on general 
2D polygonal domains. 

The linear Stokes equations are first considered in the case d = 2, ft = (0, 1) X 
(0, 1), v = 1, and / is taken to satisfy (|3.ip with a solution equal to 

(1) (^(1)^(2)) = -9(2)^(1), x (2)) 

#(i(V 2) ) ^d^^ix^Kx^) 

p(xW,xW) =100 ({x^f + {x&f) , 

denoting by ^>(x^\x^') = 1000 [x^(l - x^xj 2 '{l - x^)} 2 . The approximate 
solution (u,p) is computed with the scheme (|3.7p . The observed numerical order of 
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convergence, considering the norms ||u — P-d^\\l 2 (U) an d \\p — PvP\\l 2 (ci), is equal to 
2 for the velocity components, and to 1 for the pressure in the cases of non uniform 
rectangular and square meshes (from 400 to 6400 grid blocks). Note that in these 
cases, there is apparently no need for a significant positive value of the stabilization 
coefficient A. The observed numerical order of convergence is similar in the case of 
triangular meshes (from 1400 to 5600 grid blocks), but values such as A = 10 -4 , a = 1 
have to be used in order to avoid oscillations in the pressure field. This confirms that in 
the case of triangles, the approximate pressure space is too large to avoid stabilization. 
In fact, other tests were performed (e.g. the classical backward step) which show that 
stabilization is also needed in the case of rectangles when more severe problems are 
considered. Note that in industrial implementations, stabilization may be performed 
with other means, see [28), p], (see also [6] in the triangular case). 

We then proceed to a similar comparison in the case of transient nonlinear prob- 
lems. Considering a transient adaptation of the above steady-state analytical solu- 
tion, the continuous problem is then defined by zero initial and boundary conditions, 
T = 0.1, and the function / is taken to satisfy (|1.1[) with a solution equal to 

uW(xW,x<- 2 \t) =-t d^(x^,x^) 
u&(xW,xW,t) =tdW*(xW,xW) 
p{x^,x^,t) = 100 t ((x«) 2 + (x< 2 )) 2 ) , 

with the same function \I r as above. We again observe an order 2 of convergence of 
the approximate solution at times t = .05 and t = .1, when the space and the time 
discretizations are simultaneously modified with the same ratio (from St = 0.01 to 
St = 0.0025 as the size of the mesh is divided by 4). Similar observations arc still valid 
for the classical Green- Taylor example. 

6. Conclusions. The above numerical results show that the theoretical error 
estimate which is proved in Section [3] for the linear Stokes equations is non optimal; 
a sharper estimate is currently being written [20] under more regularity assumptions 
on the mesh. 

The proof of convergence of the full space-time discrete approximation of ()l.ll) 
given by (|4.14p uses estimates on the time translates, which were introduced in the 
L 2 (p, x (0, T)) framework for the proof of convergence of the finite volume method for 
degenerate parabolic equations [19, 15 and used for several other cases, see e.g.|2lj. A 
major difficulty which arises here is the handling on the nonlinear convective term, as 
in the continuous case, which leads us to establish an estimate on the time translates 
in L 1 (0, T; L 2 (fl)). This new technique may be used for parabolic problems with other 
type of nonlinearities. 

We remarked that industrial codes use other types of stabilizations than the 
one used here. Further works will be devoted to the mathematical study of such 
stabilizations, for which, to our knowledge, no proof of convergence is known up to 
now. 

Finally, let us also mention undergoing work on a generalization of the scheme 
studied here to the full transient Navier-Stokes equations including the energy balance, 
under the Boussinesq approximation. 
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